QCD meson spectrum in the large Nc limit 
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The low energy mass spectrum of QCD in the large Nc limit is computed. The low-lying states in the meson 
spectrum are explicitly evaluated up to the fourth order in the strong coupling perturbative expansion. The 't 
Hooft limit is smooth and the meson masses are in very good agreement with the experimental values. 



1. INTRODUCTION 

Since the seminal work of 't Hooft |1I2| the 
large Nc limit has played an increasingly impor- 
tant role in studying gauge theories in the contin- 
uum, on the lattice (see for example |3| and refer- 
ences therein) and the duality between gauge and 
string theory However, there are still no esti- 
mates of physically observable quantities, where 
the large Nc limit can be tested against experi- 
ments. Here we present the results of an explicit 
and successful calculation of the meson masses, 
in the large Nc limit In jS] we extend to an 
arbitrary number of colors the procedure used by 
Banks et al. (2) in their, by now famous, evalu- 
ation of the meson spectrum and then take the 
large Nc limit. The 't Hooft limit is smooth [7j 
and all the divergent terms in the thermodynamic 
limit cancel against the ground state vacuum en- 
ergy. 

The theory is studied in the strong coupling 
limit using staggered fermions and the Hamilto- 
nian approach to lattice gauge theory. We use 
x = l/g 2 as the expansion parameter and the 
continuum limit is extrapolated using Pade ap- 
proximants. The computation is performed by 
taking as the ground state that of the antifer- 
romagnetic Ising model which is gauge invariant 
and breaks chiral symmetry. The gauge invariant 
eigenstates of the unperturbed Hamiltonian are 
used as a basis for the perturbative expansion. 
A similar approach has been successfully used to 
study the meson spectrum of some simpler toy 
models like the Schwinger models |S] and the two 



dimensional 't Hooft model 9 . 

2. HAMILTONIAN FORMULATION 

We consider the Hamiltonian formulation of 
lattice gauge theory where the time is continuous 
and the space is discrctized on a 3-dimensional 
cubic lattice with sites x — [x\, x%, X3), where Xi 
are integers. The Hamiltonian reads 

H = Ya V £ E<1[XM2 + 7 £ Tr{UUUU) 

\ [x,n] [x.ij] 

+ J2v(n)^ A (x + n)U A B[x,n}^ B (x) + h.c. (1) 

[x.n] j 

{TrUUUU = Tr(Yl plaq[xJtj] U [xA )) where 

V(x) = (-l) z r,{y) = {-iy v(z) = (-l) v 

are the Dirac matrices for staggered fermions 
and a is the lattice spacing. The gauge field 
U[x, n] is associated with the link [x, n] and it is a 
group element in the fundamental representation 
of SU(N C ). The electric field operator E a [x, n] 
is defined on the link and it obeys the algebra 

[E a [x,n],E b [y,j]] = if abc E°[x,n]5([x,n] - [y,j]) 

and E[x,n] = E a [x,n]T a , with T a , a = 
1, N c — 1 the generators of the Lie algebra of 
U(Nc)- It generates the left- action of the Lie al- 
gebra on U[x, n] 

[E a [x,n],U[y,j]] = -T a U[x,n]S([x,n] - [y,j]) 
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The fermion fields W are defined on the lattice 
sites and obey the relation 

{^ A (x),^ B (y)} = 6 AB S(x - y) 

Each term in the Hamiltonian is gauge invariant 
under static gauge transformation. The Hamil- 
tonian is also invariant under some discrete sym- 
metries which are used to construct the opera- 
tors generating the mesons from the vacuum jl 1| . 
They are: lattice translation by a single link, shift 
along a face diagonal, parity, G-parity and Charge 
conjugation. In particular, the invariance under 
translation by a single link plays the role of a dis- 
crete chiral invariance of the theory. This sym- 
metry is broken by an explicit mass term in the 
Hamiltonian. 

2.1. The strong-coupling expansion 

In the strong-coupling expansion H is treated 
as unperturbed Hamiltonian, while Hi and H 2 
are the perturbations. Hi is the interaction 
Hamiltonian between the quarks and the gauge 
fields and Hi is the magnetic Hamiltonian: 

2a * — ' 

[x,n] 

Hi = — J2 Ti{n)^ A {x + n)U A B\x,n]^B{x) 

[x,n] 
+ h.c. 

H 2 = ttV Y TrOJUUU) + h.c. 
2g 2 a 

The vacuum of the theory is a singlet of the 
electric field algebra. It is needed then to find 
the gauge invariant eigenstates of Hq with the 
lowest energy E n . The lowest eigenstates of H 
with Eq = are singlets of the electric field al- 
gebra and since they are gauge invariant they 
must be color singlets. These states must also 
conserve charge, that means J2 X P( X )\® >= 
where p is the fermion number operator p{x) = 
*t(x)*(x) - N c /2. There are M\/{M/2)\ de- 
generate ground states. The degeneracy can be 
resolved by diagonalizing the first nontrivial order 
in perturbation theory, that is the second 



where <,>= JJ[ X n ] J dU[x,n](, ) is the inner 
product in the full Hilbert space of the model. 
This energy has been computed by constructing 
an eigenstates of Hq and using it to evaluate the 
functions present in the strong coupling expan- 
sion. If a state > is a singlet of the electric 
field algebra E a [x,n}\^ >= then H \^ >= 
and 

Q 2 

H U[x,n]\* >= Z-C 2 (Nc)U[x,n]\y > (3) 
Za 

where C 2 is the Casimir operator of SU(Nc). Af- 
ter integration over the link variable U, the vac- 
uum energy becomes 

4 2) = -K < 0\J2[p(x+n)+^][-p(x)+^](A) 

[x,n] 

where K = 1/ g 2 aC 2 Nc- The Hamiltonian has 
the two degenerate ground states of the antifer- 
romagnetic Ising model; choosing one of them 
breaks the chiral symmetry ^21- The vacuum en- 
ergy has been computed up to the fourth order in 
the strong-coupling expansion 

3. THE MESON SPECTRUM 

In the strong coupling expansion the lowest- 
lying states in the meson spectrum are those con- 
sisting of a quark and antiquark at opposite ends 
of a single link. For a given meson the wave- 
function may be determined by taking the quark 
bilinear with the desired transformation proper- 
ties, writing it in point-separated lattice form us- 
ing the discrete symmetries of the theory and ap- 
plying it to the vacuum. The mesons we have 
considered are ttq, p, u>, bi, ai, f 2 and /o- All 
the meson are degenerate at the lowest order and 
their energy is E$ — g 2 C 2 /2a. The meson en- 
ergy has been computed up to the fourth order 
in the perturbative expansion using a procedure 
similar to the one used in the evaluation of the 
vacuum energy. If one computes the meson en- 
ergy up to the fourth order, one gets 



E M 



g 2 C 2 

2a 



E 



E 



(4) 



E {2) =<0|fl a |0> + <0|.H 1 . 



n„ 



H, 



■Hi\0> (2) 



(where a is the lattice spacing). If one rescales 
the coupling constant according to the 't Hooft 
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prescription g 2 Nc — > g 2 (large A^ with g 2 Ac 
fixed) one gets for the meson masses: 

f r l 



[-+6y-171y 2 ] 



E,,, — En 



a^y 4 

' [- + 6y - 171 y 2 } 



a^/y 4 

vrip = Ep — Eq = —^[7 + 6y - 203y 2 



m bl — -E^ — E = 



l cn — E ai — Eq 



1 4 



1 A 



m fo = E fo - E 

where y = 1/g 4 . 
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1 4 



1 4 



V^ 4 



lOy - 267y 2 ] 
- 14y - 1435y 2 ] 
Uy - 875y 2 } 
18y - 1083y 2 ] 



3.1. Lattice versus continuum 

The series are valid only for g 2 large. To com- 
pare the results of the strong-coupling expansion 
with the continuum theory we need to extrapolate 
these series to the region in which y » 1. This 
region corresponds to the continuum limit be- 
cause <? 4 a 4 — > when y — > 00. One may consider 
the mass ratios, expand them as power series in y 
and then use [1, 1] Pade approximants by writing 
the mass ratios in the form P 1 = (1 + ay) / (1 + by) 
where a and b are determined by expanding to 
order y 2 and equating coefficients. In the con- 
tinuum limit this ratio yields a/b. The results 
obtained with this method are 
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-> 0.97 



For each mass ratio considered here, the [1, 1] 
Pade approximant exists with positive values for 
a/b. Therefore, an extrapolation from y = to 
y — > 00 is singularity free in this approximation. 
The results are in very good agreement with ex- 
periment except for the pion mass: this should 
not surprise due to the lack of a full chiral sym- 
metry in the lattice theory. 
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